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Using field-angle, temperature, and back-gate-voltage dependence of the weak anti-localization 
(WAL) and universal conductance fluctuations of thin Bi1.5Sbo.5Te1.7Se1. 3 topological-insulator sin- 
gle crystals, in combination with gate-tuned Hall resistivity measurements, we reliably separated 
the surface conduction of the topological nature from both the bulk conduction and topologically 
trivial surface conduction. We minimized the bulk conduction in the crystals and back-gate tuned 
the Fermi level to the topological bottom-surface band while keeping the top surface insensitive to 
back-gating with the optimal crystal thickness of ~100 nm. We argue that the WAL effect occurring 
by the coherent diffusive motion of carriers in relatively low magnetic fields is more essential than 
other transport tools such as the Shubnikov-de Hass oscillations for confirming the conduction by 
the topologically protected surface state. Our approach provides a highly coherent picture of the 
surface transport properties of TIs and a reliable means of investigating the fundamental topological 
nature of surface conduction and possible quantum-device applications related to momentum-locked 
spin polarization in surface states. 
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I. INTRODUCTION 

Similar to an ordinary-band insulator, a topological 
insulator (TI) has a bulk energy gap in its band struc- 
ture, which is generated by a strong spin-orbit interac- 
tion. The topological phase transition, brought about 
by the band inversion in the material, induces Dirac- 
fermionic surface-conducting channels!^! This topologi- 
cally protected surface state (TSS) has a helical spin tex- 
ture that is robust to small perturbations conserving the 
time-reversal symmetry, and thus prohibiting backseat - 
tering by nonmagnetic impurities!^" 8 -! 

Diverse transport studies were conducted to character- 
ize the TSS. In general, however, as-grown TIs are n- or 
p-doped so that the surface conduction can be predom- 
inated by bulk conduct ion EKD Efforts have been made 
to reduce the bulk conduction by tuning the Fermi level 
(Ep) into the bulk band gapP^H^Even w j^- n these efforts, 
however, critical inconsistencies were present in the pre- 
vious transport measurements. For instance, in bulk TIs 
with a thickness larger than ~ /im, two-dimensional (2D) 
Shubnikov-de Haas oscillations (SdHO) were observed. 
Nonetheless, the weak anti-localization (WAL) effect, rel- 
evant to the TSS, was often absent in the corresponding 
measurements, or, if present, did not fit well to the 2 D 
Hikami-Larkin-Nagaoka (HLN) WAL expressio n 1 15 1 23 ^ 
(see Appendix [A]) . In a TI, the WAL effect is generated 
by a strong spin-orbit interaction and the consequent de- 
structive interference between two electron waves trav- 
eling along a diffusive closed path in a time-reversal 
manner E2HM1 

These inconsistencies between the 2D SdHO and WAL 
were also ob served in t hinner flakes, with a thickness less 
than ~ ^ 1 21 1 22 1 35 1 j sl Furt hermore, previous 2D-SdHO 
observation j 15 | 22 | 25|)28 1 35 1 36 1 may not have b een full y rele- 
vant to the surface conduction by the TSS) 27 1 37 ^ Accu- 



rately identifying the Berry-phase shift associated with 
the TSS requires measurements in very strong magnetic 
fields, with careful Landau-level indexing! 3 ^ In most of 
the previous studies, however, the 1/2 Berry-phase shift 
was determine d based on obse rvations in relatively weak 
magnetic fieldsP^ 12 - 2 1 26 1 28 1 35 1 3 ^ Ambipolar characters with 
back-gating were also observed in the transport of TIs, 
which were assumed to be associated with the TSS. 
Here, however, the WAL effect was absent in the sam- 
ples with relatively high carrier densities ^he WAL 
effect observed in some of these ambipolar-transport sam- 
ples were reported to arise from the coupling between 
the surface a nd th e bulk bands, rather than the TSS 
exclusivelyPHHIaiHm 



It is an extremely difficult task to reliably separate 
the TSS from other conductance contributions. In this 
study, we minimized the bulk conduction using high- 
quality Bi 1 . 5 Sbo.sTe 1 . 7 Sei.3 (BSTS) TI single crystals, 
with Ep lying in the bulk gap without gating. We con- 
firmed that the WAL effect and universal conductance 
fluctuations (UCF) indeed arose from the top and bot- 
tom surfaces. By back-gate tuning the WAL character- 
istics, we identified the TSS conducting characteristics 
and the coupling between the TSS and the topologically 
trivial two-dimensional electron gas (2DEG) states that 
emerged due to band bending near the bottom surface. 
The ambipolar Hall resistivity of the bottom surface was 
consistent with the back-gate-voltage (Vb g ) dependence 
of the longitudinal resistance of the TSS. This study pro- 
vides a reliable means of differentiating the TSS of TIs 
from those of the bulk conducting state and the topolog- 
ically trivial 2DEG states, along with a highly coherent 
picture of the topological surface transport properties of 
TIs. 



2 



II. SAMPLE PREPARATION AND 
MEASUREMENTS 

BSTS single crystals were grown using the self-flux 
method 1121121 Stoichiometric mixture of high-purity start- 
ing materials (Bi(5N), Sb(5N), Te(5N), Se(5N)) were 
loaded in an evacuated quartz ampoule, which was then 
heated up to 850 °C. After annealing at 850 °C for 2 days 
to enhance the material homogeneity, the melt mixture 
was slowly cooled down to 600 °C for a week. Before 
complete furnace cooling it was kept at 600 °C for one 
more week to further improve the crystallinity. The stoi- 
chiometry and the high crystallinity of the single crystals 
were confirmed by the energy dispersive spectroscopy and 
the x-ray diffraction, respectively. 

The bulk transport properties were examined using 
^100 ^m-thick cleaved bulk crystals. For detailed char- 
acterization of transport properties with back-gating, 
BSTS flakes, which are 22 to 230 nm in their thickness, 
were mechanically exfoliated onto a Si substrate capped 
with a 300-nm-thick oxidized layer. This was then fol- 
lowed by standard electron (e)-beam patterning and e- 
gun evaporation of Ti/Au (10 nm/100— 350 nm thick) 
bilayer electrodes and contact leads. For thick crystals, 
the electrode contacts were prepared using silver paste. 
In total, four thick bulk crystals and six thin flakes were 
investigated using standard lock-in measurements, vary- 
ing T from 290 to 4.2 K. 



III. RESULTS AND DISCUSSION 

A. Thickness and temperature dependence of 
resistance 

The T dependence of the resistivity p xx of the thick 
bulk crystals of BSTS in Fig. [lja) exhibits conventional 
semiconducting behavior down to ~ 40 K. A fit of p xx (T) 
to the Arrhenius law renders the activation energy of 
E a = 26.1, 21.3, 31.6 and 20.7 meV for samples Bl, B2, 
B3 and B4 (inset in Fig. [lja) corresponds to sample B3), 
consistent with previous studies!^ However, the resis- 
tance is saturated for T below ~ 40 K, which indicates 
the emergence of additional conducting channels. This 
behavior was more pronounced in the thin flakes. Figure 
|TJb) shows a clear semiconductor-metal transition as the 
thickness of the flakes decreases. The variation of p xx {T) 
with the flake thickness can be interpreted in terms of 
surface-conducting channels in the presence of a bulk in- 
sulating gap, as illustrated in Fig. TJf). With E F in- 
side the bulk energy gap, the residual bulk conduction by 
carriers thermally activated from an impurity band was 
dominant in the thick crystals (Fig. [lja)). Thin flakes, 
however, with less bulk conductance, exhibited metallic 
behavior. One can confirm this behavior by modelling 
the simple form for total sheet conductance as follows: 
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FIG. 1. (Color online) (a) T dependence of the bulk crystal 
resistivity of the TI. Inset: Arrhenius-law fitting for the sam- 
ple B3. (b) T dependence of the normalized resistance of thin 
flakes, (c) Thickness dependence of the sheet conductance of 
thin flakes at 290 K (triangles, red online) and 4.2 K (squares, 
blue online). Solid lines are best fits to Eq. (JTJ. (d) Schematic 
measurement configuration and (e) the optical image of the 
sample F4 with invasive lead contacts. The scale bar is 2 
pm. (f) Schematic band structure near a surface of our TI 
samples. The crossed lines (red online) represent the TSS. t a 
represents the range of surface band bending (or the range of 
carrier accumulation) on the surface. Horizontal thick lines 
at t a (blue online) represent the 2DEG formed at the surface 
due to the surface band bending. Double-parabolic curves 
(blue online) are Rashba-split bands of the 2DEG. Horizontal 
dashed line corresponds to the Fermi level. 



where G s is the surface sheet conductance, Ub is the bulk 
conductivity, and t is the thickness of crystals. Here, G s 
includes the conduction through the 2DEG layer (see Fig. 
[TJf ) ) in the potential well formed by surface band bend- 
ing, as well as the conduction by the TSSpSES Fitting 
the observed results to Eq. fill (Fig. [lie)), Ob is esti- 
mated to be 86.9 and 26.7 (eyh) /znr^at 290 and 4.2 
K, respectively. These values are at least two orders of 
magnitude smaller than the ones reported previously for 
Bi 2 Se 3 p23 indicating that our BSTS single crystals were 
highly "bulk-insulating" . Assuming the range of surface 
band bending at the surface to be ~30 nm (see Appendix 
|B| in sample F4, the relative weight of the bulk to the 
surface conductance becomes <Jbt/G s ~ 26% (6%) at 290 
K (4.2 K). 
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FIG. 2. (Color online) (a) Sheet resistance Rq vs B and (b) 

Ra vs Bcos6 for different field angles at V bg = -20 V for the WAL effect Applying V bg , we confirmed that the 



served only i n the low-field range of B within a fraction 
of tesla.™ 

The 2D nature was identified more quantitatively from 
the T dependence of the MR. Figure [2jci) is the T de- 
pendence of WAL effects and the best fits of AG XX (B) 
to Eq. ([2]), from which we obtained the T dependence 
of the phase relaxation length as shown in Fig. [2jf ) 
(more details of the WAL effect are discussed below). 
Figure [2je) shows the T dependence of SG, with the cor- 
responding T dependence of the UCF amplitude <5G rms 
shown in the inset of Fig. |2jf). In a 2D system with a 
sample dimension of L ^> 1$, 1$ scales as T~ 05 for inelas- 
tic scattering by electron-electron interaction, and <5G rms 
is proportional to l^W^hi Fig.[2^f), both 1$ and SG lms 
scale as T" 5 , in good agreement with the theoretical 
predictions, indicating that the dominant inelastic scat- 
tering in the surface-conducting channels of our BSTS 
flakes was due to the electron-electron interaction. 



C. Back-gate dependence of WAL 

Up to this point, results from our BSTS consistently 
indicate that the bulk conduction was negligible, and that 
both WAL and UCF had a 2D nature. The WAL in the 
TSS arose from the Berry phase tt caused by the heli- 
cal spin texture. Since the Rashba-split 2DEG has the 
momentum- locked spin helicity (see Figs. (3^d), (e), and 
(f)), the topologically trivial 2DEG states also exhibit 



sample F4. (c) SG = G(B) - G (B) vs Bcos8 extracted 
from the data in (b), where Go(B) is the background of G(B) 
(—l/Ra(B)). For clarity, each curve is shifted vertically by 
0.03 e 2 /h. Inset: Schematic measurement configuration, (d) 
AG XX (B) = G XX (B) - G XX {B = 0) vs B at different T for 
sample F3. Solid lines are best fits to Eq. 121. (e) T depen- 
dence of SG of sample F3 extracted from the data in (d). For 
clarity, each curve is shifted vertically by 0.05 e 2 /h. (f) T de- 
pendence of the phase-relaxation length 1$ obtained from (d). 
Inset: T dependence of the root-mean-square of SG extracted 
from the data in (e). 



B. Angle and temperature dependence of WAL 
and UCF 



The surface-dominant conduction at low T becomes 
more evident in the field-angle dependence of the magne- 
toresistance (MR). Figure [2jb) shows that all of the MR 
curves taken at different field angles (Fig. |2ja)), plot- 
ted as a function of the normal component of the field 
(BjJ, merge into a single universal curve (see Appendix 
[C]for the discussion on the MR feature in in-plane fields; 
6>=90°). Even the positions of the UCF peaks agree with 
each other when plotted as a function of B x (Fig. gc)). 
These features strongly indicate that the MR in our sam- 
ple was almost completely dominated by surface conduc- 
tion over the entire field range of our measurements. Pre- 
viously, the cos(#) angle dependence of the MR was ob- 



bgi 

WAL effect arose from surface conduction, in both TSS 
and the topologically trivial 2DEG, with negligible bulk 
conduction. According to the HLN theory, for a 2D sys- 
tem in the symplectic limit, i.e., in the limit of strong 
spin-orbit coupling (r^ ^> r so ,r e ; is the dephasing 
time, t so the spin-orbit scattering time, and r e the elas- 
tic scattering time) with a negligible Zeeman term, the 
magnetoconductance correction is given as follows: 




(2) 

where ip is the digamma function, e is the electronic 
charge, h is Planck's constant divided by 2ir, and 1$ is the 
phase-relaxation length.^! Because the WAL effect con- 
stitutes a prominent transport property of the TSS, the 
relationship between the parameter a and the number of 
conducting channels in the symplectic limit is essential 
to differentiating the transport nature of TIsP^ Each 2D 
conducting channel in the symplectic limit contributes 
0.5 to the value of a. If there are two independent 2D 
conducting channels in the symplectic limit, a = a% + 
(a,, corresponding to the channel i) and is replaced 
by the effective phase relaxation length (see Appendix [D| 
for details of the WAL fitting). 

We confirmed that the back-gating affected only the 
bottom-surface conductance for the 85~90 nm-thick 
samples (F3 and F4) (see Appendix [e|. Figure^a) shows 
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FIG. 3. (Color online) (a) AG XX vs B in color codes (online) as a function Vb g . (b) AG XX (B) curves for different values of 
Vb g . Solid lines are best fits to Eq. B. For clarity, each curve is shifted vertically by 0.6 e 2 /h. (c) Schematic band diagram 
relevant to the thin flakes used in this study. The crossed lines (red online) represent the TSS. ii° p (ta 0t ) represents the range of 
the carrier accumulation on the top (bottom) surface. Horizontal thick lines at i*° p and t„ ot (blue online) represent the 2DEG 
formed at the top and bottom surface due to the surface band bending. Double-parabolic curves (blue online) are Rashba-split 
bands of the 2DEG. Horizontal dashed line depicts the Fermi level. Regions I, II, and III represent the band structure for the 
corresponding regions denoted in (g). (d, e, f) Schematic diagram of the Rashba-split strength of 2DEG for different band 
bending. Up and down arrows indicate the spin texture, (g) Vbg dependence of a, obtained from best fits to Eq. Q. Solid 
curves are the Vbg dependence of Rq. Inset: Vbg dependence of l^, also obtained from best fits to Eq. Q2p . 



AG XX vs B in color codes (online) as a function Vbg- 
Here, the WAL effect occurs over the entire range of Vb g 
of this study with a maximum AG XX at Vb g ~ —19 V, 
the Dirac point of the TSS at the bottom surface (cor- 
responding to the center diagram in Region II of Fig. 
|3jc)). Figure [3jb) shows AG XX curves for different val- 
ues of Vbg, which agree well with Eq. ^ (solid curves) 
over the entire range of B; the corresponding values of 



a are plotted in Fig. [3 
indicating that more t 



For all Vb g , a exceeds unity, 
ihan two 2D conducting channels 
with the symplectic-limit behavior were involved in the 
surface conduction. 

In Region II of Fig. |3jg), the TSS in the bottom sur- 
face contributes a value of 0.5 to a. This leaves a ~ 1 
for the top surface, which does not appear to be affected 
by Vbg. Thus, we infer that the band bending near the 
top surface is like what is shown in Fig. [3|c). In the 
top surface, in addition to the TSS, the two Rashba-split 
channels in the trivial 2DEG layer also exhibit WAL in 
the symplectic limit BUBS However, the magnitude of a 
is reduced from 1.5 (=0.5x3) to ~1 due to inter-band 
scattering, where the degree of reduction depends on the 
scattering strength.^ In Region I, Ep also enters the bulk 
conduction band (BCB) of the bottom surface. But, if 
the surface band bending is not enough to make a suf- 
ficient Rashba splitting in the 2DEG states as in Fig. 



|3jf), the band structure of the 2DEG would be simi- 
lar to the unitary casej^ where the scattering between 
the TSS and the topologically trivial 2DEG states is en- 
hanced along with weakening of the WAL effect f^D This 
reduces the value of a of the bottom surface down to 
~ 0.2 — 0.3, while leaving a unchanged at ~ 1 for the top 
surface. If Ep is shifted deeper into the conduction band 
as to form a 2DEG on the bottom surface with a large- 
Rashba-split bulk subband (Fig. ^d)), the WAL effect 
will be enhanced again, with the value of a larger than 
0.5 as shown in Fig. |3jc) for the top surface.^ In Region 
III, a similar reduction of a is expected for the bottom 
surface, due to the enhanced scattering between the TSS 
and the bulk valence band (BVB). Thus, the variation of 
a with Vbg in Fig. [3^g) is the result of variation of the 
WAL in the bottom surface state. 

The WAL e ffects report ed previousl y on TIs with 

o. ~ o.^msmm or a „ 1 mmnmMM con tained a a- 

nite bulk contribution, a ~ 0.5 corresponded to an 
effective single layer formed by the bulk and the two 
(top and bottom) surfaces, which are strongly coupled 
together. Meanwhile, a ~ 1 corresponded to an effective 
single layer formed by the n-type bulk strongly coupled 
to the top surface, in association with the p-type bottom 
surface that was decoupled from the bulk by the forma- 
tion of the depletion layer for a large negative value of 
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Vb g - To the best of our knowledge, no previous reports 
have shown good fits to the symplectic-limit expression 
of Eq. pi for fields u p to several tesla, with a exceeding 
unity \TMEmSmMSm Although the good fits of our results 
to Eq. ([2| without the Zeeman correction may be related 
to the recent report of small Lande g factor in Tlsj 38 * 56 ^ 
more studies are required to draw a definite conclusion 
on the issue. 



D. Back-gate dependence of Hall resistivity 

From the thickness, field-angle, and temperature de- 
pendence of the resistance, we conclude that in our BSTS 
samples the electronic transport was dominated by the 
top and bottom surfaces. In this case, the Hall resistivity 
can be described by a standard two-band model aa^ 
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Here, rii and Mi are the density and mobility of the car- 
riers, respectively, in the z-th conducting channel. The 
top (i=l) and bottom (i=2) surfaces constitute parallel 
conducting channels, with m being positive (negative) 
for n-type (p-type) carriers.^ In sufficiently strong fields, 
R xy converges to R xy = R% lone B sa - e ( ni +n 2 ) • In weak 
fields (B w T), Eq. ^ is reduced to 
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Here, the double signs are of the same order. The upper 
(lower) sign corresponds to ni,n 2 > (n 2 < < n{). 

Figure |4^a) shows the results of analysis of the Vb g 
dependence of Hall resistivity from sample F4. For Vt g > 
10 V, the difference between the square (red online) and 
the circle (blue online), which means nonlinearity of R xy , 
is very small, thus the R xy curves are almost linear in B. 
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FIG. 4. (Color online) (a) Vbg dependence of the slope of 
the tangent to R xy (B) curves for Sample F4 at zero magnetic 
field (circles, blue online) and 6 T (squares, red online), re- 
spectively. The curves in right inset (orange online) and left 
inset (cyan online) represent the characteristic schematic fea- 
ture of nonlinear R xy curve when the sample is in the n-n state 
As the V bg is lowered to a negative value, the R xy curves for V bg > -19 V and in the n-p state for V bg < -19 V. The 

vertical dashed line represents the boundary between the n-n 
and n-p states, (b, c) Vbg dependencies of the carrier density 
and the mobility in the bottom and top surfaces of sample 
F4, which are obtained by fitting the R xy data to Eq. pj. 
The dashed lines in (b) correspond to \n\ = 10 13 cm -2 . The 
inset in (b) shows the representative Hall resistivity curves 
for Vbg = 30 (linear dependence), —10 (n-n state, circle), and 
—30 V (n-p state, square). The solid lines are the best fits to 
Eq. {3}. 



starts to bend and the nonlinearity of the R xy 
as the Vbg decreases (curve in right inset in Fig. j4^a) ) . In 
this region (e.g. Vb g = —10 V), the slope of the tangent to 
R xy at B = T is larger than that at B = 6 T. However, 
the feature is reversed for Vb g = —30 V (curve in left 
inset in Fig. |4^a)) and the nonlinearity of R xy decreases 
as the Vb g decreases. Using Eq. Q, it turned out that the 
change in the shape of the R xy curve (from right inset to 
left inset in Fig. Qa)) indicates the ambipolar transport 
of Dirac fermions between the n-n state (top: n-doped, 
bottom: n-doped) and the n-p state on TI surfaces. 

If ^1 « [i 2 and ni,n 2 > in Eq. Q, then i?™ cak 
R s H IODS , corresponding to the region of Vb g >15 V in Fig. 
|4]ja). Since the R xy curves are almost linear in J5, the 
carrier mobility is estimated to be /xi, fx 2 ~ 140 cm 2 /(Vs) 
using the relationship /x = — , which agrees with previous 
reports. ^ 3 * 28 * 38 ! In this region (Region I in Fig. [3](c) ) with 



Ep in the BCB, the mobility decreased due to the en- 
hanced inter-band scattering.^ As Vb g decreased, with 
E F shifted to the TSS in Region II in Fig. ^c), the 
mobility of the bottom surface was enhanced so that 
Mi M2- In this case, if ni,n 2 > 0, Eq. Q leads to 
which corresponds to the curve for 
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Vbg = —10 V in the inset of Fig. |4^b) (right inset in Fig. 
Qa)). Decreasing Vbg further, Ep shifted to a p- type re- 
gion at the bottom surface. With n 2 < and \ni\ ^> \n 2 \, 
Eq. (jl) leads to |i?^ cak | < |i?^ rong |, corresponding to the 
curveTor Vb g — —30 V in the inset of Fig.Qb) (left inset 
in Fig. [4^a.) ) . The change in the relative magnitude of 
the slopes of the tangent to R xy (B), i.e., |i?^ rong | and 
|-fiW- eak |, f Qr y b ^ cross i n g —19 V clearly indicates ambipo- 
lar transport of the Dirac fermions between the n-n and 
n-p states on the TI surface. For Vb a < — 50 V (Region 
III in Fig.[3]jc)), with Ep in the BVB, scattering between 
the TSS and the BVB was enhanced once againP^ The 
resulting suppression of /i2, combined with an increase of 
ri2 in the range of Vb g < — 50 V along with the relation- 
ship (i2 = ri2efj,2 for the bottom surface, may explain the 
low sensitivity of i?n to Vb g in Fig. [3](g) . 

Fitting the R xy data to Eq. ^ gives more quantitative 
support for the analysis above on the Vb g dependence of 
the Hall resistivity. The inset of Fig. Qb) shows the rep- 
resentative Hall resistivity for Vb g = 30, —10, and —30 V, 
where the solid curves are the best fits to Eq. ^ with the 
parameter values summarized in Figs.[4^b) and (c). n\ in 
Fig.Qb) is almost constant for all values of Vf, g , while ri2 
changes its sign between n and p types at Vb g ~ — 19 V. 
This indicates ambipolar transport for the bottom sur- 
face with varying Vbg across the Dirac point, while the 
top surface remained mostly unaffected by back-gating, 
consistent with earlier qualitative analysis of Vbg depen- 
dent Hall resistivity. This back-gating effect on the two 
surfaces was also confirmed by the mobility change. In 
Fig. Qc), the best-fit values of [i\ are almost insensitive 
to the variation of Vbg. However, fi2 turns out to be sig- 
nificantly larger than n\ in the region, —50 V< Vb g <0 V, 
where Ep is assumed to be in the Dirac band of the bot- 
tom surface. The p,2 enhancement possibly stems from 
the mobility increase as Ep shifts into the Dirac band 
of the bottom surface from the trivial 2DEG band (ei- 
ther conduction or valence), where /12 is reduced by the 
scattering between the TSS and the trivial 2DEG band. 
It should be noted that, with the invasive configuration 
of electrodes adopted in this study, the observed Hall 
voltage is bound to be underestimated. However, the 
qualitative Vbg dependence of the parameters in Eq. ^ 
remains valid. 



IV. CONCLUSION 

The 1/2 Berry-phase shift in SdHO is often adopted 
to examine the topological nature of surface transport. 
However, very strong magnetic fields of B > 50 T with 
careful Landau-level indexing, required for accurate de- 
termination of the Berry phase, have made it difficult 
to clearly differentiate the conductance by the TSS from 
that by the trivial 2D-conducting states. Observation of 
SdHO also requires relatively high mobility with a suf- 
ficiently long mean-free path to support the cyclotron 
orbital motion. In contrast, the observation of WAL, an 



intrinsic 2D effect, directly points to conduction by the 
TSS. Furthermore, WAL, which arises from the coher- 
ent diffusive motion of carriers, is not limited to the high 
mobility state. In this sense, the WAL effect which was 
used primarily in this study can be considered to be a 
more essential criterion than the SdHO for confirming 
the conduction by the TSS. 

For flakes significantly thicker than an optimum thick- 
ness of ~80— 90 nm, the bulk conductance cannot be ne- 
glected. On the other hand, as the range of band bend- 
ing near the top and bottom surfaces begins to overlap 
for thinner flakes, independent gate control of the sur- 
face conduction would no longer be possible. Thus, our 
approach of separating the TSS by examining the trans- 
port characteristics specific to the 2D-topological nature 
in the optimal-thickness crystal flakes (in combination 
with back-gating) provides a convenient means of inves- 
tigating the fundamental topological nature of the sur- 
face conduction and the quantum-device applications as- 
sociated with momentum-locked spin polarization in the 
surface state of TIs. 
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Appendix A: Possible formation of multiple parallel 
2D conducting channels in TIs 

The weak anti-localization (WAL) in bulk topological 
insulator (TI) single crystals and thin T I flakes with high 
carrier density was reported previously! 22 * 25 * 31 -^ However, 
the magnitude of the consequent conductance correction 
(AG) was larger than our results by one or two orders of 
magnitude. Since the magnitude of AG is proportional 
to the parameter a in Eq. ^ in the main text, which 
corresponds to the number of parallel conducting chan- 
nels, one may suspect that multiple two-dimensional (2D) 
conducting channels connected in parallel were present 
for the conduction of TI in previous studies. A recent 
reporfP^l supports the inference. In Ref. [37] , it was con- 
cluded that the observed quantized Hall effect and SdHO 
were not caused by the topologically protected surface 
state (TSS) but by many topologically trivial 2D con- 
ducting channels connected in parallel. 

From the SdHO measurements, one can obtain the in- 
formation on the dimensionality and carrier density of the 
conducting channels. In the SdHO analysis, the degener- 
acy "2" corresponds to the bulk band or the topologically 
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trivial two-dimensional electron gas (2DEG) on the sur- 
face accumulation layer, while the degeneracy ^"cor- 
responds to the TSS. In some previous studies j^ESl the 
carrier density was estimated from the SdHO data adopt- 
ing the degeneracy '1' under the assumption that the 
observed SdHO arose from the TSS. The carrier den- 
sity estimated in this way was claimed to be relevant 
to the TSS, based on the fact that, with Ep lying in 
the TSS, the maximum carrier density is expected to be 
0.5 ~ 0.8 x 10 13 cm' 2 depending on the TI materials 
used. However, if the SdHO had arisen from the topo- 
logically trivial 2D conducting channels, the degeneracy 
should have been '2' with a doubled carrier density. In 
this case, however, a Dirac cone cannot accommodate all 
the carrier states estimated with the degeneracy '2' in 
Ref. [57] and [SU] without the bulk conduction band or 
2DEG states. 

In fact, the SdHO frequencies themselves obtained in 
Ref. [37] and Refs. (27j [60] were not much different from 
each other. Thus, the difference in the carrier densi- 
ties between Ref. [37J and Refs. [27J [SU] resulted from 
the different degeneracy values adopted in the analysis. 
Depending on the degeneracy value used in the SdHO 
analysis, one may reach very different conclusions on the 
topological nature of the conducting channels involved in 
the SdHO data. In this sense, observation of the SdHO 
itself cannot confirm the existence of the TSS. Correctly 
identifying the Berry phase in strong magnetic fields is 
essential to confirming the TSS in TIsPO 



Appendix B: Surface band bending 

The surface band bending effect is a common feature 
of semiconductors. In particular, for narrow-gap semi- 
conductors, the transport and electronic contact prop- 
erties are strongly affected by the surface band bend- 
ing. The materials which are identified as TIs are, in 
general, narrow-gap semiconductors whose band gap is 
about 100 - 300 meVP Since the energy levels of the sur- 
face state can be shifted up to a few hundred meVj^HMH 
the surface band bending has a large influence on trans- 
port properties of TIs. But, it has not been studied in 
depth to date. 

The depth of the surface accumulation layer (t a in Fig. 
[TJf ) in main text) depends on the distribution of the local 
carrier density along the z-axis! 62 * 63 ! For samples with 
the relatively high carrier density, i.e., if Ep lies in the 
bulk c onduction band, t a was calculated to be ~ 10 — 
25 nm! 20 * 22 * 35 * 45 * ^ t a c an increase further as the carrier 
density decreases! 62 * 63 * Since, in our sample, Ep lies in 
the bulk band gap with a low bulk carrier density, t a can 
be longer than 25nm. 

In addition, in comparison with the bottom surface, 
the top surface is more exposed to chemicals and e- 
beam irradiation through the sample preparation pro- 
cesses. From the careful analysis provided in the main 
text, we concluded that these processes caused the band 
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FIG. 5. (Color online) (a) The MR at V bg = -20 V for dif- 
ferent angles 9 = 0° (black online) and 9 = 90° (blue online) 
of sample F4. (b) An expanded view of the MR at $ = 90° . 



bending at the top surface, which was larger than that 
at the bottom surface, as illustrated in Fig. [3^c) in the 
main text. 



Appendix C: In-plane field dependence of 
magnetoresistance 

Figure [5] shows the MR at 9 = 90° ; direction of mag- 
netic field is in parallel with the top and bottom surfaces 
of the sample and perpendicular to the current direction 
(see Fig. [2|a) in main text). If the conduction in our 
thin flakes was only through the two surfaces (top and 
bottom) and only the localization effect affected the MR, 
the MR should have vanished at 9 = 90°. As shown in 
Fig. [5j however, a small but finite MR exists at 9 = 90° . 

The simplest inference is that the MR at 9 = 90° is 
a bulk component. In Ref. [35], the MR proportional to 
~ B 2 at 9 = 90° was observed. In the data analysis, 
this component was subtracted from the MR obtained in 
other field angles. Since the samples used in Ref. [IS] 
had a large carrier density, the large weight of the bulk 
conductance was reasonable with the ~ B 2 classical be- 
havior of the MR supporting that analysis. 

However, in our samples, as shown in Fig. [5] we did 
not find a valid argument to consider the MR at 9 = 90° 
as the three-dimensional (3D) bulk contribution. The 
~ £? 2 -type classical MR was absent at 9 = 90°. Instead, 
the MR behavior was reminiscent of the WAL effect. But, 
there is no consensus yet on whether the magnetocon- 
ductance (MC) correction (AG) of bulk carriers in TIs 
should fo llow the WAL or the weak localization (WL) 
behavior! 53 * 54 ! Thus, it is not clear whether the WAL-like 
AG(9 = 90°) in our data is of bulk origin. 

If the MR at 9 = 90° corresponds to the 3D bulk con- 
tribution, in order to extract AG of the surface conduct- 
ing channels, one has to use AG (6 = 0°) - AG(9 = 90°) 
rather than AG(9 = 0°) as used in the main text. But 
the bulk origin of AG (9 = 90°) is not clear. On the 
other hand, the magnitude of AG at 6 = 90° is suffi- 
ciently smaller than that at 9 = 0° so that the discussion 
on the angle dependence of MR and the Vb g dependence 
of MR at 9 = 0° in the main text is not affected even 
without subtracting AG(9 = 90°). Therefore, we used 
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FIG. 6. (Color online) Sets of graphs of (i) the 
digamma function — — tp I | + \ „ ) , (ii) the HLN function 



In 



function — In 



1 , ft 

2 + i^~B 



and (iii) the logarithmic 



of Eq. (SI) for i = 1, with a = 0.5 



and 1$ = 100 nm ~ 200 nm. 



the raw data for analysis of the gate dependence of WAL 
effects in Fig. [3] in main text. 

It is not clear what caused this finite MR at 9 = 90° 
in our TI flakes. It may have arisen from the side- wall 
surfaces of the thin crystal or even the in-plane MR of 
the surface conducting channels. There are some theo- 
retical prediction of in-plane field-dependence MC correc- 
tion for a 2D system, but not in the symplectic caseP^^ 
To the best of our knowledge, however, in-plane field- 
dependence MC correction of 2D systems in the sym- 
plectic limit has not been studied yet. 



However, if Z^i ^ l<j>,2, the number of fitting parameters 
increases up to 4 with a larger standard error. We solved 
this problem by taking the following simple approxima- 
tion. 

In Fig. [6j the set of curves (i) represents the digamma 
function part, the set (iii) corresponds to the logarithmic 
function part, and the set (ii) corresponds to the sum of 
the two parts. Each function is plotted with a = 0.5 and 
1$ = 100 ~ 200 nm. As displayed in Fig. [6j the digamma- 
function part is almost constant except in the weak-field 
region for different values of 1$. Thus, the HLN expres- 
sion is mostly determined by the logarithmic part. The 
digamma function causes a constant shift of the logarith- 
mic function and removes the logarithmic divergence in 
zero field. Based on this fact, four parameters in Eq. 



(Dl I can be reduced to two parameters as follows. 



Let's define l$ t i is the phase relaxation length of chan- 
nel i (i — 1,2) with the corresponding coefficient a.i 
and Z^ ff is the effective phase relaxation length with 

min-f^i,^} < Vf < max{Z^ ) i,Z^ ) 2}. Applying the ap- 
proximated behavior of the digamma function leads to 



AG. 
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Appendix D: Weak anti- localization analysis 

Magnetoconductance (MC) correction of a 2D system 
in a symplectic limit can be expressed as Eq. ([2| in the 
main text (HLN function). If there are two independent 
2D conducting channels, the equation is expanded as fol- 
lows: 



AG XX = 



i=l,2 



2n 2 h 



In 




where tp is the digamma function, e is the electron charge, 
oii corresponds to the channel i with the p hase relaxation 
length Z^J^Ilf l^x — l^_ 2 = If, Eq. (Dl I is simplified as 
follows: 



AG XX = (ai + a 2 ) 



2ir 2 h 



In 



46/2^ 



1 



4el 2 B j 
(D2) 



2ir 2 h 



ol\ In 



a 2 In 



Ael 2 a B 



(D3) 



The logarithmic part in Eq. (D3) becomes 

h 
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FIG. 7. (Color online) Three curves (blue, black and red 
online) correspond to graphs of Eqs. (D6), (D7), and (D8) 
for Qi = 1, <xi = 0.5, /^i = 100 nm and l$ t 2 — 200 nm. Inset 
shows the expanded view in the low-field range. 



where K ff 



Therefore, with Eq. (D3|, 



the Eq. (Dl ) can be simplified as 



AG X 



Q 



2ir 2 h 



In 



4e{lf) 2 B 



1> 



1 h \ 

2 + Ae(lf) 2 B J 

(D5) 



respectively. 
ln(Z^i) = In 



Here, 



4eii S 



iKW = ^(i + 4^) 



As displayed in Fig. 



the de- 



with a — ai + a.2- Figure [7] shows the validity of this 
approximation. 

In Fig. [7j the three curves (blue, black and red online) 
correspond to the followings 

aiH^x) - V(^ 2 )] + a 2 [ln(Z 0)2 ) - ^(^,2)] (D6) 
«ipn(i*,i) " ^(^,1)] + "aKW - ^(^,2)] (D7) 
aipn^,i) - + a a pn(^ l2 ) - V>(^,i)L (D8) 



and 



viation caused by different Ia in digamma function can 
be recognized only in low fields. Furthermore, since 



has a value between 1^,% and 1^,2, the deviation may be 
smaller than differences displayed in Fig. [7j Therefore, 
even with 4 parameters in different two channels, we can 
apply the one-channel HLN function with two parame- 
ters and the determined a and Is can be understood as 
a = a\ + 0L2 and Ij, = V? as Eq. 



(D5) 



Appendix E: Vbg independence of the top-surface 
conductance 

Figure [8] shows the resistance variation of an 87-nm- 
thick BSTS flake (thickness of this flake is almost iden- 
tical to that of the samples F3 and F4) as func tions of 




FIG. 8. (Color online) Vb g and Vt g dependencies of the resis- 
tance for the 87 nm-thick BSTS flake. Between two adjacent 
curves Vb g is varied by 5 V. 



back-gate {Vbg) and top-gate (V tg ) voltages. This sam- 
ple is not referred to in the main text. Except for the 
parallel shift in the resistance, the Vt g dependence of the 
resistance curves in Fig.[H]remains unaltered with varying 
Vbg. Even the positions of the resistance spikes arising 
from the UCF effect do not change for different values of 
Vbg ■ This feature indicates that the top-surface (bottom- 
surface) conductance is almost completely independent 
of Vb g (V tg ). Since this flake and the samples F3 and F4 
are of almost identical thickness we expect that the top- 
surface conductance of the two samples was independent 
of Vbg, the fact of which is utilized in our analysis in the 
main text. 
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